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Abstract 

A sample of about 5500 — > D*^ and 700 D*^a^ events is identified, using the technique 

of partial reconstruction, among 22.7 million BB pairs collected by the BABAR experiment at the 
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1 Introduction 



The neutral B meson decay modes B D*~^h~ , where h~ is a hght hadron (tt", p~, a^), have been 
proposed [|| for use in theoretically clean measurements of the Cabibbo-Kobayashi-Maskawa Q 
unitarity triangle parameter sin(2/? + 7). Since the time-dependent CP asymmetries in these modes 
are expected to be of order 2%, large data samples and multiple decay channels are required for a 
statistically significant measurement. The low efficiency of reconstructing the produced in the 
D* decay leads to significant loss of signal events. Partial reconstruction of the B meson results in 
substantially larger efficiency, albeit with higher backgrounds. The overall sin(2/? + 7) sensitivity 
is expected to be roughly similar for partial and full reconstruction, and about 90% of the partially 
reconstructed events cannot be fully reconstructed. Therefore, both full and partial reconstruction 
can and should be used for the sin(2/3 + 7) measurement. 

The measurement of the B^ lifetime, presented here, constitutes a first step toward measuring 
sin(2/3 + 7) with partial reconstruction of B^ — > D*^ . In this analysis, we have developed the 
procedures for candidate reconstruction, background characterization, vertexing, and fitting, all 
critical components of the time-dependent sin(2/3 -|- 7) analysis. 

2 The BaBAR Detector and Data Sample 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. 
The data consist of 22.7 million BB pairs, corresponding to an integrated luminosity of 20.7 fb~^ 
recorded at the T(45) resonance. In addition, 2.6 fb~^ were collected about 40 MeV below the 
resonance. This off-resonance sample is used to study the continuum, e~'~e~ — > qq background, 
where q = {n, d, s, c}. 

The BABAR detector, described in detail elsewhere |^], consists of five sub-detectors. Surround- 
ing the beam-pipe is a five-layer silicon vertex tracker (SVT), providing precision measurements 
of the positions of charged particles close to the interaction point and tracking of charged par- 
ticles with low transverse momentum. Outside the support tube that surrounds the SVT is a 
40-layer drift chamber (DCH), filled with an 80:20 helium-isobutane gas mixture. The DCH pro- 
vides charged particle momentum measurements in a 1.5 T magnetic field, and ionization energy 
loss measurements that contribute to charged particle identification. Surrounding the DCH is a 
detector of internally reflected Cherenkov light (DIRC), providing charged particle identification. 
Outside the DIRC is a CsI(Tl) electromagnetic calorimeter (EMC), used mainly to detect and 
measure the energy of photons and to provide electron identification. The EMC is surrounded by 
a superconducting coil, which generates the magnetic field for tracking. Outside the coil, the flux 
return iron is instrumented with resistive plate chambers, used mainly for muon identification. 

PEP-II is an asymmetric energy storage ring, with positron and electron beam energies of 
about 3.11 and 9.0 GeV. The center-of-mass (CM) frame of the average e~^e~ collision is therefore 
boosted along the z direction in the lab frame, enabling time-dependent measurements involving 
reconstruction of B mesons. 

3 Analysis Method 



8 



3.1 Partial Reconstruction 



To partially reconstruct & ^ D*^p~ candidate^ , only the p and the vr^, the soft pion from the 
decay D*~^ D^'^t decay, are reconstructed. The angle between the momenta of the B and the 
p in the CM frame is then computed: 



Ml 



+ EcuEp 



cos 9 Bp 



2Pb\Pp 



(1) 



where Mx is the mass of particle x, Ep and pp are the measured CM energy and momentum of the 



p, Eqm is the total CM energy of the beams, and Pg = y Eq-^/4: — M^o- Given cos 6 Bp and the 

measured four-momenta of the tTs and the p, the B four-momentum may be calculated up to an 
unknown azimuthal angle (p around pp. For every a value of (p, the expected D four-momentum, 
VDi4'), is determined from four-momentum conservation, and the (/>-dependent "missing mass" is 
calculated, m{(p) = \/|'Pd(0)P- With mmax and m^i^ being the maximum and minimum values of 
m((/)) obtained by varying cj), we define the missing mass, mmiss = ^ [mmax + m-mm]- This variable 
peaks around Mj^o, the nominal mass, for signal events, with a spread of about 3.5 MeV, while 
background events are more broadly distributed. 

We define the D* helicity angle 9b)* to be the angle between the directions of the D and the B 
in the D* rest frame. The value of cos6d* is computed as in Ref. The p helicity angle 6p is 
defined as the angle between the directions of the vr'^ (from the decay of the p) and the CM frame 
in the p rest frame. Since the B^ D*~^p~ decay is (87.8 ± 4.5)% longitudinally polarized the 
cos Bp and cosOd* distributions of signal events peak toward ±1. Background events have fairly flat 
distributions, with BB background increasing in the region of soft vr^ and tt", making cos9d* and 
cosOp useful for background suppression. With knowledge of cos ^d* and cos 9bp, the 3-momentum 
of the D is determined, up to a two- fold ambiguity and detector resolution effects. 

3.2 Event Selection 

In order to suppress the continuum background, we select events in which the ratio of the 2nd to the 
0th Fox- Wolfram coefficients Q is smaller than 0.35. Charged p candidates are identified by their 
decay to a "hard" charged pion, vr^, and a vr*^. To suppress fake vr^ candidates, the vr*^ momentum 
in the CM frame is required to be greater than 400 MeV. The invariant mass of the vr" candidate 
must be within 20 MeV of the nominal vr^ mass. Both the vr/j and vr^ candidate tracks are required 
to originate within 1.5 cm of the interaction point in the x — y plane (the plane perpendicular to 
the beam) and within ±10 cm of the interaction point along the direction of the beams. Tracks 
are rejected if their specific ionization and/or Cherenkov angle indicate that they are highly likely 
to be a kaon or a lepton. The invariant mass m(p) of the p candidate must be between 0.45 and 



The computed cosine of the angle between the B and p must satisfy | cos Bspl < 1. To suppress 
combinatoric background, we require |cos0p| > 0.3 and |cos0£)*| > 0.3, and reject events in the 
range cos^p > 0.3 and cos6'd. < —0.3. 

For each p^Trf pair that satisfies the above requirements, m^iss is calculated. The pair with 
the smallest value of | mmiss — in the event is selected, and all others are discarded. Similarly, 

the "wrong-sign" p^vr^ pair with the smallest | mmiss — -^d"! is retained for background studies, 
as described below. Right-sign events in the range 1.810 < mmiss < 1-840 GeV are classified as 

^Charge conjugate decays are also implied. 




1.1 GeV. 
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sideband events, and are used for background studies. Right- and wrong-sign events satisfying 
?^miss > 1-845 GeV are classified as signal region events. All other events are discarded. The 
efficiency of signal events to satisfy all the signal region criteria is 6.4%. 

Events that could be fully reconstructed in the decay modes K~tt^ or K~7r~^7r^ are 

tagged as fully reconstructed. In addition to satisfying the partial reconstruction criteria above, 
fully reconstructed events are identified by requiring that the reconstructed D invariant mass be 
within 40 MeV of Mj^o , the difference between the D* and D invariant masses be between 142 and 
150 MeV, the reconstructed CM B energy be within 50 MeV of £^cm/2, and uies > 5.25 GeV, 

where tties = \/^cm/^ ~ P\ beam energy substituted mass and is the reconstructed CM 

momentum of the B meson. 

3.3 Mecisurement of the Decay Time Difference 

The decay position Zrec of the partially reconstructed B candidate along the beam direction is 
determined by constraining the tt/j track to originate from the beam-spot in the x — y plane. To 
account for the B meson flight in the x — y plane, 30 /um are added in quadrature to the beam-spot 
size. The tTs is not used in this vertex fit in order to simplify the classification of the background, 
and since its contribution to the vertex precision is small, due to multiple scattering. 

The decay position -Mother of the other B meson along the beam direction, is obtained with all 
tracks excluding the vr/j, the tTs, and any track whose CM angle with respect to cither of the two 
calculated directions of the D is smaller than 1 radian. This "cone cut" reduces A'^^, the number 
of D daughter tracks used in the other B vertex. The remaining tracks are fit with a constraint to 
the beam-spot in the x — y plane. The track with the largest contribution to the of the vertex, 
if greater than 6, is removed from the vertex, and the fit is carried out again, until no track fails 
this requirement. 

The decay time difference At = (zrec — Zother)/j(3c is then calculated, where 7/3 is the CM frame 
boost. The value of 7/? is continuously determined from the beam energies, and averages 0.55. The 
estimated error a^^ in the measurement of At is calculated from the parameters of the tracks used 
in the two vertex fits. 

Events are rejected if the probability of the tt/j vertex fit is smaller than 1%, or if the 
probability of the other B vertex is smaller than 0.5%. We also require \At\ < 15 ps and 
0.3 < cT^j < 4 ps. 

The quantities Zrec and Mother are computed in the same way for fully reconstructed candidates 
as they are for partially reconstructed candidates. Using the Monte Carlo simulation, it is verified 
that with the 1 radian cone cut, the TV^^ distribution of events that are fully reconstructed in the 
mode R- "Tr"*" or K 7T~^7T^ is in good agreement with the distribution of partially reconstructed 
events. 

3.4 Backgrounds 

The types of events in the partially reconstructed on-resonance sample are classified as follows: 

1. Signal: B^ —>■ D*~^p~ events, in which the tt/j is correctly identified. This requirement ensures 
that Zrec is the decay position of the B meson, up to the effect of detector resolution. The tTs 
or the 7r° candidates may be mis-reconstructed. 

2. B° ^ D*+ai: 5° ^ D*+a'[ events, in which the TTh is a daughter of the oi, and hence 
originates from the decay point of the B meson. 
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3. Peaking B^B^ background: B^ — > D*^p^ and some B^ — > L'*^a^ events, in which the iTh 
originates from the other B meson, resulting in the measurement At = 0, up to the effect 
of detector resolution and the selection of tracks used in the other B vertex. The mmiss 
distribution of these events peaks around Mj^o, similar to signal events. 

4. "Combinatoric" BB background: Random combinations of tt/j, , and vTs candidates, possi- 
bly including true p decays. 

5. B ^ D**p~, where D** stands for a charged or neutral resonance with mass in the range 
2.4 - 2.5 GeV decaying into D*^Tr. 

6. Continuum, e'^e^ — > qq events. 



3.5 Probability Density Function 

An unbinned maximum likelihood fit is used to obtain the B^ lifetime r „ from the data. The 
fit is performed simultaneously to on- and off-resonance data, and to on-resonance events which 
were fully reconstructed. The probability density function (PDF) is a function of At, a^^, and 
four "kinematic" variables: 1) mmiss! 2) m{p); 3) rn-Es! ^-iid 4) F, a Fisher discriminant that helps 
distinguish between BB and qq events. The value of F is computed from the total CM energy 
flow of tracks and neutral EMC clusters (excluding the p and the tTs) into nine volumes, defined by 
nine 10°-wide concentric cones centered around the p CM momentum p^^. Each cone is folded to 
combine the energy flow in both hemispheres with respect to ph- 

The PDF of partially reconstructed on-resonance events is a sum of terms corresponding to the 
different event types: 

T-'ii) = fsigna.l'PsignaliO + fD'aiVD'aiiO 
+ /peakfiO ^peakB0(O + f Bb'^ BB^O 

+ fD"VD"{i) + fm'Pm{i), (2) 

where ^ = (mmiss, ^(p), F, At, a^J is the vector of fit variables, 

ViiO = M^im^i,,) n,{m{p)) T^{F) %{At, ci^J (3) 

is the PDF corresponding to event type i, and fi is the fraction of events of type i in the data 
sample, where X]j fi = ^■ 

The PDF of the off-resonance sample is Vqq, which is also used to describe the continuum 
component of the on-resonance events in Eq. (^) . 

The PDF of fully reconstructed events is similar to V{^), except that J'i{F) is replaced by the 
function £'(mEs)- The fractions fB,*ai, fpeakB^^ f^BB ^^^^ sample are determined from the 
Monte Carlo simulation, and are of order a few percent. The fractions /fignai ^-^d /^^ are obtained 
from a 3-dimensional fit to the mmiss, m(p), and mEs distributions of this sample. 

The 

mmiss distribution of signal events is parameterized as a bifurcated Gaussian, 

( ^ ( (?"miss-M)2\ 
signal (mmiss) « exp I ^ I , (4) 

where M is the position of the peak, and the value of a depends on the sign of mmiss — The 
proportionality constant in this and subsequent PDF expressions is determined by integrating the 
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PDF over the allowed range of the PDF variable. The m^iss distributions of the background event 
types are parameterized as a bifurcated Gaussian plus an ARGUS function M , 



^("Imiss) OC mmiss\/l " (^^miss/^A) 

X exp e (l - (mmiss/MA) 



(5) 



where ^(mmiss 

from fits to data or Monte Carlo simulation, as described in Sec. 3.6 



for mjniss > Ma, and Ma and e are parameters whose values are determined 



The TZi{m{p)) functions are sums of a relativistic P-wave Breit Wigner function and second- 
order polynomials. The functions J^i{F) are bifurcated Gaussians, and £i{mEs) are a Gaussian for 
signal events and ARGUS functions for the backgrounds. 

The At PDF of signal events is an exponential decay with the lifetime, convoluted with a 
triple-Gaussian resolution function to account for finite detector resolution: 



fn Gji (tr , (7^ 



+ fw Gw(tr, <T^^) -\- fo Go{tr, (T^ 



(6) 



where At^ is the true decay time difference between the two B mesons, tr = At — Att is the At 
residual, G„, G^, and Go are the "narrow", "wide", and "outlier" Gaussians, each of the form 



1 



2'irsa. 



■ exp 



(tr- - by 



(7) 



where s and b are parameters obtained from the fit to data. The coefficients fi of Eq. (^) satisfy 
fw = ^ — fn — fo- The same At PDF parameters are used for signal and B^ D*^a^ events. 

The At PDFs of the combinatoric BB, peaking B^B^, B D**p~, and continuum are of the 
form 



rBgd.(At,cj^j 



1 



/-^e-l^**l/- + (l-/05(Att] 



1 - /o) 1 dAtt 

fn Gn{tr, CT^J + (1 — fn) Gyj{t^,a 



+ fo 



1 



2ttSo 



exp 



{U - boY 
2sl 



, (8) 



where the phenomenological parameter r is different from v^g of Eq. @. The 6{Att) term accounts 
for events in which the TTh originates from essentially the same point as the tracks dominating the 
determination of Mother • Parameter values for each of the different background PDFs are determined 
from fits to independent control samples in data, as described in Sec. The fraction fo of events 
corresponding to the outlier Gaussian is taken to be zero for the three BB background PDFs, based 
on studies performed with the Monte Carlo simulation and the data. In the continuum PDF, is 
determined to be about 1% for So = 10 ps. In the peaking B^B^ PDF fj- = 0, reflecting the fact 
that the tt^ originates from the decay of the other B. 



3.6 Fit Procedure 

The B^ lifetime t^q is obtained through a series of flts. First, the kinematic variable PDF parameters 
of B^ D*^a^ , combinatoric BB and the peaking B^B^ background are determined by fitting the 
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distributions of these variables in the Monte Carlo simulated events. The parameters of ^signaiC-^)) 
as well as fn^ai, and /peaks") a-^e also obtained from the Monte Carlo simulation. 

The Tl{m{p)), M.{mmiss)i and £{mEs) parameters of signal and continuum events, and the 
J-{F) parameters of continuum are obtained by fitting the kinematic variable distributions of the 
data in the signal region. This 4-dimensional fit is performed simultaneously for partially and fully 
reconstructed on-resonance data, and for off-resonance data. The fractions fqq and f^^ are also 
determined in this fit, with /signal obtained from /j = 1. The value of fo** is set to in this 
and the subsequent fits, and is later varied when studying systematic errors. 

The parameters of are determined by fitting the data in the rn-miss sideband. The sideband 
is populated only by combinatoric BB and continuum events. Consequently, this fit uses F as 
the only kinematic variable, and is performed simultaneously for on- and off-resonance data. The 
parameters of 7^q(At, cr^J of this sample are also determined in the fit. 

The composition of the wrong-sign event sample in the signal region and the parameters of 
T^eakfio (^^5 f'^At) determined from fits to this sample. The parameters of 7^g(At,(T^J of this 
sample are also determined in the fit. 

Finally, the signal region parameters of 7^ignai(At, cr^J and 7^q(Af,cr^J are determined from a 
simultaneous fit to the right-sign signal region on-resonance, off-resonance and fully reconstructed 
data. The parameters of Tg^( At , a^^ ) and 'Tweaks" (^^i ^i^t ) ^'^^ taken from the sideband and wrong- 
sign fits, respectively. Use of the sideband and wrong-sign samples for this purpose is validated 
using the Monte Carlo simulation. 

The samples used to obtain parameters of the different At PDF components are summarized 
in Tab. |. 

Table 1: Signal and background At PDFs and the data samples from which their parameters are 
determined. The indicated data samples are used simultaneously in the fits, where "on" and "off" 
refers to on- and off-resonance data, respectively, and "full" refers to the fully reconstructed event 
sample. 



PDF 


Data samples 


m-miss region 


Right / wrong-sign 


T^^(At,(j^J 


on + off 


Sideband 


Right-sign 


'^eakBo('^*)<7At) 


on + off 


Signal region 


Wrong-sign 


'T^ignal ( At, cr^ J 


on + off + full 


Signal region 


Right-sign 



4 Results 

The partially reconstructed signal region on-resonance sample contains 50898 events, including 
5266±251±34 W D*^p- and 691±36±4 5° D*+a^ events, as determined by the kinematic 
variable fit. The systematic errors are due to the finite numbers of events in the Monte Carlo 
simulation samples used to obtain the kinematic variable distributions of the two BB backgrounds 
and B^ — > D*~^ai . The fully reconstructed sample contains 255 it 20 signal events. 

Projections of the PDF onto the data are shown in Figs. [l| through ^. The result of the maximum 
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Figure 1: Projections of the PDF of partially reconstructed events onto the on-resonance data in 
the variables mmiss (upper left), m{p) (upper right), F (lower left), and At (lower right). Events in 
the At projection plot satisfy the additional criteria mmiss > 1-854 GeV, 0.6 < m{p) < 0.93 GeV, 
and F < -2.1. 



likelihood At fit is 



raw 



1.535 ± 0.064 (stat.) ps, 



(9) 



where the error is statistical only. 

Several corrections are applied to this result, in order to account for known sources of bias. A 
correction of +0.008ib0.011 ps is due to the assignment of oi daughter tracks to the other B vertex, 
decreasing the measured At in L)*"''a^ events. The resulting value of 1.543 ps is divided by 

Rd = 0.982 lb 0.022, the ratio between the B'^ lifetime obtained from a fit to signal Monte Carlo 
events and the value obtained from a fit to the true At^ distribution of this sample. This correction 
accounts for the effect of D daughter tracks that pass the cone cut and are used in the other B 
vertex. A correction of +0.014 it 0.013 ps is applied to this result, to account for a possible bias 
due to event selection, determined by fitting the true Atj distribution of signal events passing the 
selection criteria. The magnitudes and errors of all these corrections are obtained from the Monte 
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Figure 2: Projections of the PDF of continuum events onto the off-resonance data in the variables 
"imiss (upper left), m(p) (upper right), F (lower left), and At (lower right). 



Carlo simulation. 

The fit PDF was used to generate and fit hundreds of Monte Carlo samples, each corresponding 
to the data sample in number of events and PDF parameters. The value of obtained from these 
fits was on average lower than the generated value by 0.031 it 0.005 ps. Repeating these studies with 
different Monte Carlo sample sizes, this bias is understood to be due to limited sample statistics. 
A correction of this magnitude is therefore added to the value of t^q . The fully corrected result is 

r^o = 1.616 ± 0.064 (stat.) ps. (10) 

5 Systematic Errors and Cross Checks 

Several sources of systematic error are considered. The statistical error matrix obtained from the 
wrong-sign signal region At fit is used to vary the parameters of the peaking B^B^ background, 
taking into account their correlations. A fit to the right-sign signal region data follows each variation 
in these background parameters. The resulting variations in v^g are added in quadrature to form 
the total systematic error due to the finite number of events in the wrong-sign sample. With the 
same procedure, the errors due to the sideband fit are propagated to the wrong-sign signal region 
and then to the right-sign signal region fits, to obtain the error due to the sideband sample size. 
The errors due to the finite number of events used in the kinematic variable fits on data and the 
Monte Carlo simulation are evaluated in the same way. 

The Monte Carlo statistical errors in the determination of Rd, the B^ D*^a^ bias, and the 
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Figure 3: Projections of the PDF of fully reconstructed events onto the fully reconstructed data in 
the variables mmiss (upper left), m{p) (upper right), ttt-es (lower left), and At (lower right). 



selection bias corrections are taken into account. The fraction of events in which D daughter tracks 
are used in the other B vertex fit is varied by ±5%. The magnitude of this variation is determined 
by comparing the distributions of fully reconstructed data and Monte Carlo events. The 
resulting variation \n is used to evaluate the systematic error due to this uncertainty. 

The contribution of the B — > D**p~ background is neglected in the fits. To evaluate the 
systematic error associated with this, we instead take the number of S — > D**p^ in the data 
sample to be 2400, corresponding to B(B D** p') B{D** D*+tt) = 0.3%. This value is 
estimated from known branching fractions of the decays B — > D**Tr~ , B D^*^ p~ , B — > D^*\~ ^ 
and available limits on B{B D**p~) Q. The PDF parameters of this background are taken 
from the Monte Carlo simulation. Repeating the At fit yields a 0.023 ps change in the value of 
T^p, which is taken as the systematic error. The D** states simulated are -Di(2420), -D2(2460), and 
Di{j = ^), the latter having mass 2.461 GeV and width 290 MeV. No significant difference is found 
between these states in terms of their effect on the observable quantities of this analysis. 

Two variations of the signal At PDF are used in fits to the data. In one variation, the parameter 
b in Eq. (j^ is replaced by ba^^ . In the other, the sum of the narrow and wide Gaussians of Eq. (^) 
is replaced by a Gaussian convoluted with an exponential. The effect of generating Monte Carlo 
event samples using one PDF and fitting them using another PDF was studied. Based on these 
data and Monte Carlo studies, a systematic error of 0.016 ps is estimated due to the choice of signal 
At PDF. 

The parameters s and b (Eq. (0)) of the signal and continuum PDF outlier Gaussians are fixed 
in the fit to the right-sign signal region data. To estimate the associated systematic errors, their 
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Table 2: Systematic errors. 



Source 


Error (ps) 


Statistical error of sideband fit 


0.033 


Statistical error of kinematic fit 


0.029 


Statistical error of wrong-sign fit 


0.002 


Monte Carlo statistics: Calculation of Rd 


0.036 


Monte Carlo statistics: Kinematic parameter fits 


0.014 


Monte Carlo statistics: Event selection bias 


0.013 


Monte Carlo statistics: D*~^a^ bias 


0.011 


uncertainty 


0.026 


Level of -B ^ D** p" background 


0.023 


Likelihood fit bias 


0.016 


Variation of fixed parameters 


0.015 


B{B^ D*+a^)/B{B^ D*+ p-) 


0.005 


Level of peaking background 


0.003 


Bias from fully reconstructed events 


0.001 


SVT misalignment 


0.008 


z-length scale uncertainty 


0.007 


Beam energies uncertainty 


0.002 


Total 


0.075 



values are varied within reasonable ranges, and the resulting changes in v^g are taken as systematic 
errors. 

Additional systematic errors are due to the uncertainty in the relative branching fractions of 
B^ D*'^a^ and B^ D*^ p~ , the level of peaking background, the introduction of a possible 
bias due to the use of fully reconstructed events, detector alignment and z-length calibration, and 
beam energy uncertainty. The total systematic error is 0.075 ps, dominated by the errors due to 
sideband and kinematic fit statistical errors, and Monte Carlo statistical errors. The systematic 
errors are listed in Tab. |2|. 

Several cross-checks were conducted to ensure the validity of the result. The number of signal 
events detected is in good agreement with the published branching fraction |8| and our signal 
reconstruction efficiency. The fit was repeated with different values of the cone cut, ranging between 
0.6 and 1.2 radians. The data were fitted in bins of the lab frame polar angle, azimuthal angle, and 
momentum of the tt/j, and in sub-samples corresponding to different SVT alignment calibrations. 
In all cases, no significant variation of the result was observed. 

6 Conclusion 

In a sample of 22.7 million BB pairs, we identified 5521 ± 252 ± 34 ^ D*+ p- and 691 ± 36 ± 4 
^0 ^ events using partial and full B reconstruction. These events were used to measure 
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the lifetime, with the prehminary result being 

To = 1-616 ± 0.064 (stat.) ± 0.075 (syst.) ps. (11) 

This result is in good agreement with earlier published measurements Q, constituting a necessary 
step in validating the use of partially reconstructed D*^p~ events for the measurement of 

sin(2/3 + 7). 

Acknowledgements 

We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of this 
project also relies critically on the expertise and dedication of the computing organizations that 
support BABAR. The collaborating institutions wish to thank SLAC for its support and the kind 
hospitality extended to them. This work is supported by the US Department of Energy and National 
Science Foundation, the Natural Sciences and Engineering Research Council (Canada), Institute 
of High Energy Physics (China), the Commissariat a I'Energie Atomique and Institut National 
de Physique Nucleaire et de Physique des Particules (France), the Bundesministerium fiir Bildung 
und Forschung (Germany), the Istituto Nazionale di Fisica Nucleare (Italy), the Research Council 
of Norway, the Ministry of Science and Technology of the Russian Federation, and the Particle 
Physics and Astronomy Research Council (United Kingdom). Individuals have received support 
from the A. P. Sloan Foundation, the Research Corporation, and the Alexander von Humboldt 
Foundation. 

References 

[1] P.F. Harrison and H.R. Quinn (ed.), BABM Physics Book, Chap. 7.6 (1998); R.G. Sachs, 
Enrico Fermi Institute Report, EFI-85-22 (1985) (unpublished); I. Dunietz and R.G. Sachs, 
Phys. Rev. D37, 3186 (1988) [E: Phys. Rev. D39, 3515 (1989)]; I. Dunietz, Phys. Lett. B427, 
179 (1998). 

[2] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. Theoret. 
Phys. 49, 652 (1973). 

[3] The BABAR Collaboration, B. Aubert et al, Nucl. Instr. and Methods A479, 1 (2002). 

[4] The CLEO Collaboration, G. Brandenburg et al, Phys. Rev. Lett. 80, 2762 (1998). 

[5] The CLEO Collaboration, ICHEP98 852, CLEO CONE 98-23 (1998). 

[6] G. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[7] The ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B254, 288 (1991). 

[8] See, for example. The Particle Data Group, C. Caso et al., Review of Particle Physics, Eur. 
Phys. J. C15, 606 (2000). 



18 



